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Finite-temperature quantum criticality in a complex-parameter plane 

C. Li and Z. Son£] 

School of Physics, Nankai University, Tianjin 300071, China 

A conventional quantum phase transition (QPT) occurs not only at zero temperature, but also ex¬ 
hibits finite-temperature quantum criticality. Motivated by the discovery of the pseudo-Hermiticity 
of non-Hermitian systems, we explore the finite-temperature quantum criticality in a non-Hermitian 
PT-symmetric Ising model. We present the complete set of exact eigenstates of the non-Hermitian 
Hamiltonian, based on which the mixed-state fidelity in the context of biorthogonal bases is calcu¬ 
lated. Analytical and numerical results show that the fidelity approach to finite-temperature QPT 
can be extended to the non-Hermitian Ising model. This paves the way for experimental detection 
of quantum criticality in a complex-parameter plane. 

PACS numbers: 11.30.Er, 64.70.Tg, 05.70.Jk 


I. INTRODUCTION 

A central task of the theory of quantum phase transi¬ 
tions (QPTs) is to describe the consequences of this zero 
temperature singularity on physical properties at finite 
temperature since all experiments can only be done at 
non-zero temperature [I|. It turns out that the mixed- 
state fidelity, which is related to the statistical distance 
between two density operators, is a powerful way to de¬ 
scribe the signatures of QPTs at nonzero temperature 
0 . 

Recently, there have been intense efforts to establish a 
“PT-symmetric quantum theory as a complex extension 
of the conventional quantum mechanics , since the 

discovery that a non-Hermitian Hamiltonian having si¬ 
multaneous parity-time (VT) symmetry has a real spec¬ 
trum Q. Motivated by the pseudo-Hermiticity of non- 
hermitian systems, the QPT in a non-Hermitian VT- 
symmetric Ising model, which is driven by a complex 
staggered transverse field, has been explored based on 
the exact solution [12] . It has been shown that the phase 
diagram can be characterized by the geometric phase. 
The phase boundary can be identified by a divergence of 
Berry curvature density, which is defined in the context 
of biorthogonal bases. 

In this paper, we address two questions concerning the 
quantum criticality for the non-hermitian Ising model at 
finite temperature. At first, whether the QPT driven by 
complex parameters has residual criticality at finite tem¬ 
perature. Second, as a metric-based approach, whether 
the mixed-state fidelity in the context of biorthogonal 
bases is available to describe this criticality. Analytical 
and numerical results indicate that the phase boundary 
at finite temperature can be identified by the mixed-state 
fidelity. 

The paper is organized as follows: In Sec. [TT]we present 
the exact solution of the non-Hermitian Ising model. In 
Sec. HD the mixed-state fidelity of thermal states is ob¬ 
tained analytically and a numerical simulation is per- 


* songtc@nankai.edu.cn 


s 

^ 'At 


f 0i) 

/ ( p\ 

-l \ / 

° 1 1 


J 7 

\S r 


Ar _ 

—At 

Oil) 


FIG. 1. (Color online) Phase diagram for the ground state 
of the Ising ring in a complex staggered transverse field. The 
heavy lines (green) represent the boundary which separates 
three quantum phases. Phases I and III are paramagnet, while 
II is ferromagnet. Vector denotes an arbitrary position in 
the £ r] plane, while "r^ =L A“r^ are two points in the vicinity 
of point The mixed-state fidelity for those two points can 
characterize the thermal criticality at point "r^. 

formed to demonstrate the signatures of QPTs at finite 
temperature. Section [TV| contains the conclusion. 


II. MODEL AND SOLUTION 

In the last few years, there have been many studies of 
the phase diagrams for various quantum models in the aid 
of the fidelity j~L3l [Toj or mixed-state fidelity [!, l20l - l22| . 
All the investigations are performed for Hermitian Hamil¬ 
tonians, where the probability is preserving in the context 
of standard Dirac inner product. It would be interesting 
to explore the extension of this metric-based approach to 
a non-Hermitian Hamiltonian with pseudo-Hermiticity, 
where biorthonormal inner product is introduced. We 
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begin our investigation by introducing a non-Hermitian 
Ising model, whose ground state phase diagram has been 
studied in previous work 0. It has been shown that the 
Berry phase in the context of biorthonormal inner prod¬ 
uct can identify the phase diagram, which sheds some 
light on the extension of the metric-based approach to a 
non-Hermitian system. 

The concept is studied exemplarily for criticality in 
a non-Hermitian one-dimensional spin-1/2 Ising model 
with a complex staggered transverse magnetic field on a 
27V-site lattice. The system is modeled by the following 
Hamiltonian 


2N 

n =-JXh°i°i+i (!) 

3 = 1 

where complex field g. = rj + i (— 1)° £ (i = V^l) with 
r] and £ being real numbers. Here dj (A = x, z) are 
the Pauli operators on site j, and satisfy the periodic 
boundary condition = cr^ + 2N . In ground state this 
model exhibits quantum critical points at rf 1 + £ 2 = 1 
and r] = 0 for rj 2 + £ 2 > 1, separating paramagnetic and 
ferromagnetic phases 0- 

Following the derivation in Appendix, all the eigen¬ 
states of H can be constructed as the form 

i*>= n ( 2 ) 

kE[ 0,7r] 

with eigen energy 

E= £ e n (k). (3) 

{l”>fe} 


Here | n) k (n G [1,16]) denotes the 16-dimensional qudit 
(16-level system) state, which is expressed explicitly as 


n G [1,5] : | n) k 

= ^r[ 2cos ( fc / 2 ) 

—2i sin (k/2) 


, PkP-k 

e n + i2t e" - i2t 


1 


. a kPk a -kP-k + 


e n + 2r, kHk ~ kH ~ k ' e n -2 V 

+e ik/2 4^_ k - e- ik /^W_ k ] |0), 


( 4 ) 


and 


|6>1 = 71 ( e “ /! d/sh + e-“ ,2 d«h) l»>. 

F)* = 4/sllo>, M = °Vh|0), (5) 


and 

n G [9,12] : | n) k 

= + i(e n ) 2 -(v + i0 2 - 1] 

-2(7?cos 2 | - it sin 2 ^)}e~ lk/2 al 
+ [(e" + r?) 2 +e 2 -l]cos^ 

-i[(e n - it) 2 -rf - 1] sin ^of_ k a\0[ 

-i (t? + smke~ lk/2 alpl k pl} |0), (6) 

and 

TIG [12,16] : R, 

= ^T {{(e n + V-i£) I(n 2 -(v + i0 2 - 1] 

-2(7? cos 2 ^ - it sin 2 ^ )}e~ lk,2 a ] _ k 

+ [( e ” + v) 2 + C 2 — 1] COS 7}P-k 

-i[(e n - it) 2 ~rf - 1] sin ^ot\of_ k ^_ k 

-i (t? + it) sin ke~ %kl2 |0), (7) 

where Q n is corresponding normalization factor. And 
coefficients e n (k) are given as 

e 1 = [2r 2 cos (2ip) + 2 + 2 \/r A — 2r 2 cos k + l] 1 / 2 , 

e 3 = [2r 2 cos (2 tp) + 2 — 2\/r 4 — 2r 2 cos k + l] 1 / 2 , 
e 2 = —e 1 , e 4 = -e 3 , e 5 = e 6 = e 7 = e 8 = 0, 
e 9 = [r 2 cos (2(p) + 1 

+ \/4r 2 cos 2 ip + 2r 2 cos k — 2r 4 ] x / 2 , 
e 11 = [r 2 cos (2y?) + 1 

-74^2 cos 2 + 2r 2 cos k — 2r 4 ] 1 / 2 , 
e 14 = -e 13 = £ 10 = -6 9 , 

c 16 = -el 5 = C 12 = -6 n . (8) 

Here a\ and are fermionic creation operators and we 
parameterize the complex field in terms of the polar ra¬ 
dius and angle 


r = \/r\ 2 + £ 2 and tan cp = %/rj. (9) 

Similarly, we can obtain the eigenstates of 1-0 by the 
qudit state (n| fe , which has the form 


MOU = (l n (-?))ft) t - ( 10 ) 

It is turned out that the complete biorthonormal set can 
be constructed due to the fact 

( n \k I m )k' = 5 mn&kk', I n )k H k = 1 - ( n ) 
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FIG. 2 . (Color online) 3-D plots of mixed-state fidelity of finite-temperature thermal state of the non-Hermitian Ising model. 
We consider the case with N = 300 and the field variances Ar] = A£ = 0.01J -1 . Figures (a-c) plots of F as functions of 77 
and £ for different temperatures /3, while (d-f) are corresponding inversions of (a-c). At high temperatures (a), the mixed-state 
fidelity shows a flat surface. When the temperature decreases (b, c), the decay of fidelity becomes sharp and clearly shows the 
signature of QPT which occurs at absolute zero temperature. A more accurate dependence is plotted in Fig.3. 


At zeroTemperature, the phase diagram is obtained 
in Ref. 0 and sketched in Fig HI The influence of 
the imaginary field on the phase transition is obvi¬ 
ous: It enhances the action of the real transverse field 77 , 
shrinking the disorder region, leading to a circular phase 
boundary. In the following, we will show that this circle 
critical line is still robust at finite temperature. 


III. MIXED-STATE FIDELITY 

The key point of this paper is that the mixed-state 
fidelity 


F (p, p) = tT^/p 1 / 2 pp 1 / 2 ( 12 ) 

can be utilized to characterize the thermal criticality al¬ 
though density matrices p and p are non-Hermitian. The 
non-Hermiticity of the matrices arises from the descrip¬ 
tion of two different thermal states 


p = exp (—/3T-L) /tr exp (—/3T-L ), (13) 

p = exp (~PH) / tr exp (-PH) , (14) 

where 

H = H(r l ,0,H = U(r 1 + 5r l ,t + Sr 1 ), (15) 

are both non-Hermitian operators. 

As the consequence of the translational symmetry, 
[hi c, h k '] = 0 , leads to 


F (p, P) = n tr \/ d/VdA ( 16 ) 

k 

which can simplify the computation procedure. Studies 
on Hermitian systems show that a zero-temperature fi¬ 
delity vanishes at the quantum boundary. It is expected 
that the similar behavior can occur for the present non- 
Hermitian model. Actually, in the limit case f3 oo, 
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F (p, p) is reduced to 


Foo(p,p)=n o k=i[\(i\ k \i) k 

k k 


(17) 


by using the completeness condition in Eq. CU). Here 
states (l\ k and \l) k are 


(4 = <1 T + A?)| fc , |l) fc = |1(V-AV)>„ (18) 

where denotes a vector in the plane (see FigITj) . 
In the following, we consider a simple case with T* = 
(cos p, sin ip) and AT" = (Ar cos p, Ar ship). Corre¬ 
spondingly, states (l\ k and \l) k can be expressed in a 
polar coordinate system (r, <p) 


<l| fc = (1 (1 + Ar, <p)\ k , |l) fc = |1 (1 - Ar, <p)) h . (19) 

We note that any vanishing overlaps (l| fc |l) fe can result in 
the zero point of fidelity F^ (p, p). In the thermodynamic 
limit N — * oo, k becomes a continuous variable. This 
ensures us to estimate the overlap lim^o (l| fe |l) fe . Direct 
derivation shows that 

lim |1 (r, <p)) k = - fiW-k 

+ 7 * a W-k + 7 PlP-k) | 0 ), (r < 1 ), ( 20 ) 

and 

lim (l(r,<p)\ k ~ (0|i, (r > 1), ( 21 ) 

k — 

where (p E [0, U (^,27r], and 7 = ^/l — r 2 sin (p + 
ir sin (p. Meanwhile we have 


for Ar > 0, which is reduced to yjAr/2 for Ar <C 1. In 
this sense, points (l,<po) are a little special compared to 
other points on the circle. This may be due to the fact 
that points (l,po) are three-phase points. In this paper, 
we focus on the critical points on the circle except these 
two points. 

For finite temperature, numerical simulations are per¬ 
formed for different /3. We plot F as a function of (£, rj) 
in Fig 12] and Figj3l We see that F is almost a constant 
and equals to unity for a wide range of (£, 77 ), apart from 
the very narrow area around the circle p 2 + £ 2 = 1, where 
drastically drops to zero as shown in Figj3l For larger /?, 
we can see that the dip of F well indicates the critical 
point of the non-Hermitian Ising model for different val¬ 
ues of £. We also study the temperature dependence of 
the fidelity for different critical fields, which is labeled by 
(p. In Fig. 0]we plot F as functions of temperature for 
several typical values of (p and N = 10 4 . We fit the data 
with the function 


In F = T/3 + In A, (26) 

which indicates the fidelity decays exponentially with de¬ 
cay constant T and the pre-exponential factor A. Here T 
and A are functions of (p as the form 

T = [-1.1 + 1.5 cos (2<p)] x 10 5 , 

In A — [—7.1 + 9.5 cos 2(p — 2.4 cos (4<p)] x 10 11 . (27) 

Based on above analytical and numerical results, we con¬ 
clude that the metric-based approach established for Her- 
mitian systems can be applicable to the non-Hermitian 
Ising model. 


lim |1 (r, <p)) k = -4_(a41 ft - p\a ] _ k 

+7 + 7 PlPU) 1 °). (r < 1 ), ( 22 ) 

and 

\im {l(r,cp)\ k = {0\ia k /3 k a-k/3-k, (r > 1), (23) 

k—^0 

where tpG (f ? ^ L )- These lead to the conclusion that 


lim (1 (1 + Ar, ip) U1 (1 - Ar, <p)) fc = 0, (24) 


for Ar > 0 and (p ^ ^, ■^ L . The essence of the results 
is the nonanalytical behavior of states lim/^o (l| fe and 
lim/c^o |l) fe at the critical points. This indicates that 
the fidelity approach for QPT can be extended into the 
complex regime at zero temperature. 

However, the overlap in both points <po = \ and ^ 
can be calculated analytically as 


lim 

k-rO 


(1 (1 + Ar, ¥> 0 ) I*, I 1 (! - Ar, <p 0 )) k 
Ar 


(1 + Ar) y^2Ar — (Ar) 2 


(25) 


IV. SUMMARY AND DISCUSSION 

In summary, we have shown that the QPT in a non- 
Hermitian Ising model has residual criticality at finite 
temperature, which can be described by the mixed-state 
fidelity in the context of biorthogonal bases. The success¬ 
ful extension of mixed-state fidelity to complex regime 
may be ascribed to the following two facts: At first, we 
happen to have an exact solvable non-Hermitian model, 
whose ground state is divided into several regions, or 
phases in the a complex-parameter plane. So far, this 
model is unique, while the QPT in most discovered non- 
Hermitian models is referred as the exceptional point. 
Secondly, the biorthogonal complete set is employed to 
formulate the density matrix of a thermal state. In this 
situation, we cannot give a general conclusion for other 
non-Hermitian systems, which should be investigated in 
the future. 

Our study shows that quantum criticality in such a 
non-Hermitian system appears remarkably even at finite 
temperature. This paves the way for experimental detec¬ 
tion of quantum criticality in complex-parameter plane. 
In experiment, there are several ways to realize the Her- 
mitian Ising model, such as using trapped ions p3l - [26j , an 
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FIG. 3. (Color online) 2-D plots of mixed-state fidelity of finite-temperature thermal state of the non-Hermitian Ising model. 
We consider the cases with N = 10 4 and field parameters (a) £ = 0, Ar] = A£ = 0.01J -1 , (b) £ = 0.4, Ar] = A£ = 0.01J -1 and 
(c) £ = 0.8, Ar] = A£ = 0.01J -1 . Three curves represent the fidelities at various temperatures indicated in figures, respectively. 
It shows that the decay of fidelity for low temperature well indicates the critical point of the non-Hermitian Ising model. When 
the temperature increases, the decay of fidelity becomes less sharp, leading to boundary broadening of QPT. 



FIG. 4. (Color online) Plots of F as functions of temperature 
for several typical values of Lp and N — 10 4 , where (3 is in the 
unit of . We fit the data with the functions in Eqs. ([26]) 
and (1271) . 


array of cavities mm, atoms within a cavity [29]. and 
Rydberg atoms in an optical lattice [30l-[33j. Recently, it 
was pointed that [34|, an imaginary transverse field may 
be implemented by optically pumping a qubit state into 
the auxiliary state. 


Appendix: Exact solution of the 7i 


We consider the solution of the non-Hermitian Hamil¬ 
tonian of Eq. ©• We start by taking the Jordan-Wigner 


transformation [35| 


a t =n i 1 - 2c i ci ) c p 

(A.1) 

l<3 


aj = i - 2 c]cj, 

(A.2) 

= - n i 1 - 2c i ci ) ( c i+4) ’ 

(A.3) 


l<3 


to replace the Pauli operators by the fermionic operators 
Cj. Likewise, the parity of the number of fermions 

2 N 

n = nK) = (- 1 ) Arp ( A - 4 ) 

1=1 

is a conservative quantity, i.e., [H,n] = 0, where N p = 


Y^ 2 j=i c ] c j- Then the Hamiltonian (pQ) can be rewritten 
as 


where 


U= Y, PvKaPa, (A.5) 

1(1 +an) (A.6) 


is the projector on the subspaces with even (a = +) and 
odd (a = —) N p . The Hamiltonian in each invariant 
subspace has the form 


2N—1 

H a = -J Y ( c i c i+ 1 + c ]+i c i + c i£+i + c j+i c i) 
j= i 

+Ja (c\n c 1 + 4 c 2 n + 4jv4 + c l c 2 Jv) 

2N—1 

—(l — 2c}cj) (A.7) 

j =i 

taking the Fourier transformation 


E jka ( a/c ( , even j 
... I Pk c , odd j ’ 


(A.8) 
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for the Hamiltonians 'H (T , we have 

U a = -jY,Hk„ (A.9) 

k(j 

Hk a = (e lka + l) + {e lka ~ .1) a \ a ft- kd , + H - c - 

+2r] - 2 (rj + i£) 

-2( V -iOPljK (A.10) 

where the momentum k a are defined as fc+ = 

2 (m + 1/2) 7 t/N, k- = 2m7r/N , m = 0,1, 2,TV — 1, 
respectively. 

In this paper, we focus on the case with large TV, in 
which the effect of fc+ — fc_ can be neglected. In the 
following we will neglect the subscript a in k G . In order to 
obtain the eigen states in both subspaces with even and 
odd N p , what we need is to diagonalize the Hamiltonian 

h k = H k + H_ k (A.11) 

with 

H k = 2e ik / 2 cos (fc/2) a{p k + i2e ik ' 2 sin (fc/2) a\p f _ k + H.c. 

+27? - 2 (7? + it) a{a k - 2 ( 7 ? - *0 fifa, (A. 12) 

due to the fact [h k ,h k r] oc 0. This procedure is 
similar to introduce the Bogoliubov transformation for 


solving a simple transverse field Ising model. There 
are two 8-dimensional invariant spaces for h k , spanned 
by two sets of bases {a\_ k |0), /?t k |0), a\ k aL k 0\ k |0), 

4*4*4*10)} and {|0>, 4*4*10), 4 fe /?L |0), 

44* |0), PkP-k l°)> 44 a -*4* |°)}, involving odd 
and even numbers of fermions, respectively. Here |0) is 
the vacuum of fermion operator Cj, i.e., Cj |0) = 0. For¬ 
tunately, such two 8x8 matrices can be diagonalized 
analytically in the form 

h k \n) k = 2e n (k) \n) k ,n G [1,16], (A.13) 

where \n) k and e n (k) are listed in Eqs. (4)-(10), which 
can be used to construct the eigenstates and spectrum 
of the original Hamiltonian. By the same procedure, we 
can obtain the eigenstates of TO. 
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